12.2 SCOPE vs STOP

Fit LME models for Prolactin, HPL and GH2.
Fit GEE models for Insulin and HOMA-IR.

  • Outcomes: Prolactin, hPL, GH2, Insulin, HOMA-IR, Glucose (all outcomes log-transformed)
  • Independent variables: Study, GA at sampling
  • Model fitted: LME or GEE model
          log(hormones) ~ Study + GA at sampling + Study:GA at sampling
          LME only: Random intercepts for Study ID and Plate
          GEE only: Independence correlation structure

12.2.1 Global tests

Type III ANOVA on main and interaction effects, and LM model summaries.

PRL (ng/ml)
Chisq Df Pr(>Chisq)
(Intercept) 21.782244 1 <0.0001 ***
Study 1.759986 1 0.2 ns
GA_samp 37.833421 1 <0.0001 ***
Study:GA_samp 1.825160 1 0.2 ns
hPL (μg/ml)
Chisq Df Pr(>Chisq)
(Intercept) 107.37637 1 <0.0001 ***
Study 63.37312 1 <0.0001 ***
GA_samp 115.88145 1 <0.0001 ***
Study:GA_samp 71.34390 1 <0.0001 ***
GH2 (ng/ml)
Chisq Df Pr(>Chisq)
(Intercept) 9.974225 1 0.002 **
Study 26.572090 1 <0.0001 ***
GA_samp 57.571753 1 <0.0001 ***
Study:GA_samp 26.560416 1 <0.0001 ***
Insulin (uU/ml)
model term df1 df2 F.ratio Chisq p.value
1 Study 1 Inf 3.706 3.706 0.05 ns
3 GA_samp 1 Inf 0.183 0.183 0.7 ns
2 Study:GA_samp 1 Inf 2.665 2.665 0.1 ns
HOMA-IR
model term df1 df2 F.ratio Chisq p.value
1 Study 1 Inf 4.810 4.810 0.03 *
3 GA_samp 1 Inf 0.432 0.432 0.5 ns
2 Study:GA_samp 1 Inf 2.900 2.900 0.09 ns
Glucose (mmol/L)
model term df1 df2 F.ratio Chisq p.value
1 Study 1 Inf 0.138 0.138 0.7 ns
3 GA_samp 1 Inf 0.001 0.001 1.0 ns
2 Study:GA_samp 1 Inf 0.190 0.190 0.7 ns

12.2.2 Model summary

Ratio of Geometric means (95% CI)

PRL (ng/ml)
Var Effect P
(Intercept) 2.14 (1.24, 3.03) <0.0001
StudySTOP 0.63 (-0.30, 1.56) 0.18
GA_samp 0.18 (0.12, 0.24) <0.0001
StudySTOP:GA_samp -0.04 (-0.10, 0.02) 0.18
hPL (μg/ml)
Var Effect P
(Intercept) 2.53 (2.04, 3.01) <0.0001
StudySTOP -2.02 (-2.52, -1.52) <0.0001
GA_samp 0.17 (0.14, 0.20) <0.0001
StudySTOP:GA_samp 0.14 (0.11, 0.18) <0.0001
GH2 (ng/ml)
Var Effect P
(Intercept) -1.04 (-1.70, -0.38) 0.002
StudySTOP -1.77 (-2.45, -1.08) <0.0001
GA_samp 0.16 (0.12, 0.20) <0.0001
StudySTOP:GA_samp 0.12 (0.07, 0.16) <0.0001
Insulin (uU/ml)
Var Effect P
(Intercept) 3.45 (2.16, 4.73) <0.0001
StudySTOP -1.31 (-2.64, 0.02) 0.05
GA_samp -0.05 (-0.13, 0.04) 0.27
StudySTOP:GA_samp 0.07 (-0.01, 0.16) 0.10
HOMA-IR
Var Effect P
(Intercept) 2.16 (0.83, 3.50) 0.001
StudySTOP -1.55 (-2.94, -0.16) 0.03
GA_samp -0.06 (-0.14, 0.03) 0.21
StudySTOP:GA_samp 0.08 (-0.01, 0.17) 0.09
Glucose (mmol/L)
Var Effect P
(Intercept) 1.66 (1.30, 2.02) <0.0001
StudySTOP -0.07 (-0.43, 0.29) 0.71
GA_samp 0.00 (-0.02, 0.03) 0.84
StudySTOP:GA_samp -0.01 (-0.03, 0.02) 0.66

12.2.3 Diagnostics

Model diagnositcs.

Diagnositcs for PRL (ng/ml)

Diagnositcs for hPL (μg/ml)

Diagnositcs for GH2 (ng/ml)

Diagnositcs for Insulin (uU/ml)

Diagnositcs for HOMA-IR

Diagnositcs for Glucose (mmol/L)

12.2.4 Estimated means

Estimated marginal means (or adjusted means) adjusted for GA at sampling.

PRL (ng/ml)
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 46.77984 8.557414 Inf 32.68503 66.95279
STOP 10 60.98198 1.984145 Inf 57.21453 64.99751
SCOPE 16 157.65554 4.671804 Inf 148.75981 167.08324
STOP 16 145.80292 10.047989 Inf 127.38133 166.88860
hPL (μg/ml)
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 58.66173 5.1415039 Inf 49.40261 69.65621
STOP 10 37.37513 0.7441673 Inf 35.94469 38.86251
SCOPE 16 191.79504 2.4328863 Inf 187.08546 196.62318
STOP 16 239.75276 7.7025414 Inf 225.12153 255.33490
GH2 (ng/ml)
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 1.5266388 0.1577893 Inf 1.246690 1.869452
STOP 10 0.9794978 0.0222874 Inf 0.936775 1.024169
SCOPE 16 4.6508967 0.0835049 Inf 4.490076 4.817477
STOP 16 5.2495097 0.2270129 Inf 4.822907 5.713847
Insulin (uU/L)
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 19.77659 4.6279509 Inf 12.50145 31.28545
STOP 10 11.15166 0.3842101 Inf 10.42348 11.93070
SCOPE 16 14.96218 0.6124863 Inf 13.80862 16.21210
STOP 16 13.12747 0.9877112 Inf 11.32756 15.21337
HOMA-IR
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 5.004074 1.2160703 Inf 3.107908 8.057111
STOP 10 2.360224 0.0858766 Inf 2.197770 2.534686
SCOPE 16 3.587401 0.1563751 Inf 3.293639 3.907364
STOP 16 2.735166 0.2179779 Inf 2.339631 3.197568
Glucose (mmol/L)
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 5.395297 0.3344260 Inf 4.778085 6.092238
STOP 10 4.773670 0.0366162 Inf 4.702440 4.845979
SCOPE 16 5.475695 0.0719187 Inf 5.336536 5.618483
STOP 16 4.691147 0.0561450 Inf 4.582385 4.802490

12.2.5 Comparisons

Hormone Contrast GA_samp Ratio of Geometric means (95% CI) P-value
PRL (ng/ml) STOP / SCOPE 10 1.30 (0.91, 1.88) 0.2
PRL (ng/ml) STOP / SCOPE 16 0.92 (0.80, 1.07) 0.3
hPL (μg/ml) STOP / SCOPE 10 0.64 (0.53, 0.76) <0.0001
hPL (μg/ml) STOP / SCOPE 16 1.25 (1.17, 1.34) <0.0001
GH2 (ng/ml) STOP / SCOPE 10 0.64 (0.52, 0.79) <0.0001
GH2 (ng/ml) STOP / SCOPE 16 1.13 (1.03, 1.24) 0.010
Insulin (uU/ml) STOP / SCOPE 10 0.56 (0.35, 0.90) 0.02
Insulin (uU/ml) STOP / SCOPE 16 0.88 (0.74, 1.04) 0.1
HOMA-IR STOP / SCOPE 10 0.47 (0.29, 0.76) 0.002
HOMA-IR STOP / SCOPE 16 0.76 (0.64, 0.91) 0.003
Glucose (mmol/L) STOP / SCOPE 10 0.88 (0.78, 1.00) 0.05
Glucose (mmol/L) STOP / SCOPE 16 0.86 (0.83, 0.89) <0.0001

12.2.6 Barplots by Study

Barplots in this section may be misleading due to involvement of Study by GA interaction.

12.2.7 Plots for GA